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Abstract We investigate a specific set of two-loop self-energy 
corrections involving squared decay rates and point out that 
their interpretation is highly problematic. The corrections can- 
not be interpreted as radiative energy shifts in the usual sense. 
Some of the problematic corrections find a natural interpre- 
tation as radiative nonresonant corrections to the natural line 
shape. They cannot uniquely be associated with one and only 
one atomic level. While the problematic corrections are rather 
tiny when expressed in units of frequency (a few Hertz for hy- 
drogenic P levels) and do not affect the reliability of quan- 
tum electrodynamics at the current level of experimental accu- 
racy, they may be of importance for future experiments. The 
problems are connected with the limitations of the so-called 
asymptotic-state approximation which means that atomic in- 
and out-states in the S-matrix are assumed to have an infinite 
lifetime. 
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1 Introduction 



In view of the rapid progress in ultra-accurate Lamb shift measurements in atomic hydrogen , 
it appears useful to investigate the mathematical foundations of the theorems on which our 
level-shift calculations are based, in addition to the continuing efforts (e.g. [^-Q) of evaluating 
higher-and-higher order radiative corrections to the bound-state energy levels. Here, we focus on 
a number of mathematical subtleties connected with the famous theorem of Gell-Mann, Low and 
Sucher [||,|9| that forms the basis for the derivation of the expressions investigated in level-shift 
calculations. This theorem has historically proven to be an extremely useful tool in the analysis 
of bound-state quantum electrodynamics, and it yields formal expressions for the renormalized 
level shifts whose numerical and analytic evaluation has lead to accurate predictions for the 
bound-state energy levels that are in agreement with all experiments reported so far in the 
literature. However, second thoughts about some mathematical subtleties connected with the 
theorem may be required at the level of accuracy envisaged in projected experiments. 

Consider the formula for the energy shift of an atomic state as given by the well-known theorem 
of Gell-Mann, Low and Sucher |^,[| 

i s A d 

AE n = lim ———ln[{n\S S) \\ri) c }, (1) 

where |n) is an unperturbed asymptotic bound-electron state as given by the Dirac theory. S £ \ 
is the infinitesimally damped S matrix given by the time-ordered exponential 



S e> x = Texp ^-iA J dt J d 3 x exp(— e|t|) Hi(x)^J 
where x denotes the four- vector (t, x), and the interaction Hamiltonian density reads 



(2) 
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The interaction Hamiltonian density involves the quantized electromagnetic field A^(x) , the 
quantized Dirac field ip{x) in the Furry picture, and the mass counter term 5m. The index c 
in ([l]) indicates that only connected graphs (see e.g. Ch. 6 of [pi) enter into the expression 
for the energy shift. Note that in writing down the expression (||), we implicitly assume the 
physical existence of the asymptotic state \n), i.e. of the unperturbed state \n) with an infinite 
lifetime. If the interaction with the quantized electromagnetic field (the "vacuum modes" ) could 
be "switched off" [as it is assumed for the damped interaction (||) in the distant past and 
future], then all states would be asymptotic states and could be used to construct 5-matrix 
elements rigorously. While this is - strictly speaking - unphysical, Eq. (|2|) is still an excellent 
approximation for most bound-state calculations. 

A second-order evaluation of ([!]) in powers of Ttj, which involves the one- loop self-energy, shows 
that the radiative energy shift of an excited state with nonvanishing angular momentum (say, a P 
state) has an imaginary part generated by the interaction of the atomic state with the quantized 
electromagnetic field (which gives the expression for the decay rate T). This effect limits the 
validity of the asymptotic-state approximation which is a good approximation as long as T is 
much smaller than the separation of atomic energy levels [see also the discussion in Ch. XXI 
of especially Eq. (3.28) on p. 547 ibid.). 

Certain two-loop self-energy contributions to the hydrogenic energy levels involve the square 
of the imaginary parts of two one- loop insertions, which result in a shift of the real energy 
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eigenvalue. The error initially made in ignoring the decay width of the resonances thus influences 
the real part of the eigenvalues which theory predicts at two-loop order. This leads to a problem 
in connection with the asymptotic-state approximation that is originally used in writing down 
the expression ([l]). At present, we have no better way to gauge the magnitude of this problematic 
effect but to evaluate it, according the current formalism, within a theory that a priori involves 
asymptotic states. The effects discussed here touch a certain question regarding the mathematical 
foundations of bound-state quantum electrodynamics. 

This paper is organized as follows: in Sec. [2L we present an evaluation of the problematic two-loop 
corrections involving the squared decay rates for the 2P state of atomic hydrogen. In Sec. ^ we 
argue that some of the problematic corrections find a natural interpretation as radiative correc- 
tions to the off-resonance effects that influence the line shape in atomic transitions. Conclusions 
are left to Sec. [|. 



2 A Concrete Example 



An expression for the two-loop self-energy correction to the energy of a bound hydrogenic state 
within the formalism of nonrelativistic quantum electrodynamics (NRQED) has been derived 
in 0] [see Eq. (16) ibid.], and we take this equation as the starting point of our investigations. 
Diagrammatically, the two-loop self-energy can be represented as in Fig. |l]. We will investigate 
the effect within the e-method that involves a scale-separation parameter e for the photon en- 

For the two photons with 



ergy (this method is explained in [12| and the Appendix of [0 
energies u% and L02, we need two scale-separation parameters e\ and €2- These cancel when the 
high- and the low-energy parts are added, and it is permissible to keep only the divergent and 
constant terms as both e\ — > and €2 — ► 0. We therefore keep €1/2 a s variables and evaluate 
only those contributions to the two-loop self-energy integrals that correspond to the square 
of the residues along the uj\ and u>2 integrations (these correspond to the "squared imaginary 
parts" or "squared decay rates" ) . We focus on the 2P state of atomic hydrogen and use natural 
Gaussian ("microscopic") units with h = c = 6q = 1. The Schrodinger energy of an atomic 
state is E n = —(Za) 2 m/(2n 2 ) where n is the principal quantum number, a is the fine-structure 
constant, and Z is the nuclear charge number. 

There is a contribution due to the diagram with crossed loops in Fig. |l| (a), 
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Here, the p l are the momentum operators, and H is the Schrodinger Hamiltonian. We "pick 
up" only the terms of the "squared-decay" type, i.e. the terms generated by the infinitesimal 
half-circles around the poles at uj\ = E2P — -Eis and L02 = E2P — Ei$. For the evaluation of these 
terms, the specification of the infinitesimal imaginary part — i 5 is required in order to fix the sign 
of the pole contribution. For the contribution Ci(2P) generated by the poles at u>i = E2P — -Eis 
and L02 = E2P — E\s in 7i(2P), we obtain 
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Figure 1: The crossed- loop (a), the rainbow diagram (b) with one loop 
inside the other, and the loop-after-loop diagram (c) which contribute 
to the two-loop self-energy for a bound electron. The propagator of the 
bound electron is denoted by a double line, and the arrow of time is 
from left to right. 



where the summation convention is used and the matrix element M\ reads 
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Note that the contribution C\ lacks the factors ir in the denominator which are characteristic of 
other two-loop corrections: these are compensated by additional factors of tt in the numerator 
that characterize the pole contributions. 

The rainbow diagram in Fig. p] (b) with the second loop inside the first does not create squared 
imaginary contributions. From the irreducible part of the loop-after-loop diagram in Fig. [I] (c) 
(we exclude the reference state in the intermediate electron propagator), we obtain 



T 2 (2P) 



lim 

5^0+ 



x ( 2P 



2a 
3 tt m 2 

1 



ei 



P 



H — i5 — Eiv + 



P 



(■1 



1 



H — E%p 



dW2 L02 



P> 



1 



H — i 5 — E2P + 0J2 



P' 



2P 



(8) 



Again, picking up only those terms which are generated by the infinitesimal half-circles around 
the poles at uj\ = E 2 p — Eis and u>2 = E2P — Eis, we obtain the contribution C2(2P) involving 
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squared decay rates, 
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where the matrix element M. 2 reads 
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(10) 



The prime in the reduced Green function indicates that the 2P state is excluded from the sum 
over intermediate states. 

This contribution deserves a more detailed discussion because C 2 can be written as a "second- 
order perturbation" according to 
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where Ej^(.E) denotes the one-loop nonrelativistic self-energy operator with the energy argu- 
ment E; the matrix elements of this operator read 
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In ( pT| ) and (12), we pick up only the term corresponding to the square of the two imaginary 
contributions and obtain the expression 
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where the sum covers all P states except the reference state (both discrete and continuum 
states). Now, we complete the sum over the intermediate states, observing that only P states 
yield a nonvanishing contribution. After angular averaging over atomic momenta, we obtain a 
contribution of 
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in agreement with @. 

From the derivative term (reducible part of the loop-after- loop diagram), we obtain 

2 m re 2 
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where = (2/3) 8 a (Za) 4 m is the decay width of the 2P state. The last contribution of the 
"squared-decay" type - it originates from the "seagull term" characteristic of NRQED - reads 
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Adding all contributions, we obtain a shift of 
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(19) 



for the 2P level. For atomic hydrogen {Z = 1), this correction amounts to —14.9 Hz. Compared 
to the total Lamb shift of the 2P level, for which a value of —12 638 380(80) Hz has been given 
in [|l4| (for the 2P 1/ / 2 -state), this is a tiny effect. The comparison illustrates the accuracy of the 
predictions obtained using the Gell-Mann-Low-Sucher theorem (|]) • The theoretical value of the 
2P Lamb shift has been improved in part by the recent numerical evaluation of the one-loop self 



energy for low Z [15|, but the removal of the main theoretical uncertainty will necessitate the 
complete evaluation of the two-loop self-energy corrections of order a 2 (Za) & . Currently, only 
the double logarithm is known in this order in the (Za)-expansion [|i~6| , |i~7i , and the groundwork 
for the evaluation of single-logarithmic and nonlogarithmic corrections has been laid in H]. In 



Eq. (19), we write the final result as a contribution to the i?60 coefficient (for a definition of the 
analytic ^-coefficients of the two- loop self-energy, see e.g. ]i~3"|). 

For the 3P state, we have to take into account the decays into the IS and 2S states. For example, 
the contribution Ci(3P) reads 
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where \ = (E 3 p — -E^s) (E%p — E 2 s)- The sum of C\ - C4 for the 3P state of atomic hydrogen is 
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For atomic hydrogen, this correction evaluates to —7.47 Hz. 
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3 Interpretation of the Squared Decay Rate 



In this section, we will investigate the question of whether the squared decay rates receive a 
natural interpretation within the formalism of scattering theory. 

We consider the scattering amplitude associated with the diagram in Fig. 0. The hydrogenic 
atom in the ground state is excited by a laser photon with frequency wl that is close to the 
resonance ul ~ i?2P — -Eis- Within the resonance approximation, we may restrict the sum over 
intermediate states to the 2P level only, as indicated by the label "2P" for the electron line. 
This is the so-called resonance approximation. 



cjl ^ -E-2P — -Bis 



IS S 2P IS 




(b) 



~ B 2 p - -Bis 

A 2P r s 2P r s 2P is 



Figure 2: Contributions to the photon scattering cross section near the 2P resonance 
of atomic hydrogen. The laser photon of frequency ul ~ E2P — Eis 1S absorbed, and 
the dominant contribution to the scattering amplitude is from the virtual 2P state, 
as indicated by the label on the electron line. The spontaneously emitted photon 
has a frequency u> 8 . Figs, (b) and (c) represent radiative corrections to the scattering 
amplitude. The arrow of time is from left to right. 
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Fig. g (a) represents the dominant Kramers-Heisenberg contribution without radiative cor- 
rections. In Fig. |H (b), we have a one-loop self-energy insertion in the electron propagator. We 
assume that the absorbed photon is very close to the resonance (u/l ~ E 2 p — ^as) an d set the 
energy argument of the self-energy insertion equal to E 2 p. Using these approximations, the self- 
energy insertion, within the resonance approximation, can be described by the diagonal matrix 

>L \e 2 p) 2P) of the self-energy operator (12). The imaginary part of this matrix 



element ^2P 
element reads 



where 
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is the well-known decay rate of the 2P state. The diagram in Fig. || (c) involves two one-loop 
self-energy insertions and entails in that sense the "square of the decay rate" of the 2P level 
within the resonance approximation. This square of the decay is usually not interpreted as an 
energy shift. Rather, one sums the infinite series of one-loop insertions of which the first terms 
are shown in the Feynman diagrams in Figs. | (a,b,c). We ignore in the sequel the real part of 
the self-energy operator and define the "decay rate operator" T via the relation 
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This leads to the following infinite series representing the electron propagator in the resonance 
approximation with an infinite number of one-loop self-energy insertions (of which we discard 
the real part and keep only the imaginary part — iT2p/2), 
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This result expresses the well-known fact that the decay rate term — iT2p/2 in the electron prop- 
agator denominator is generated by summing an infinite series of Feynman diagrams involving 
one-loop self-energy insertions; the first three terms in this series are shown in Figs. § (a) - (c). 

Let us now go beyond the resonance approximation inherent to Eq. (|2^ ) and consider off-resonant 
atomic states (see Fig. H). We therefore replace 



|2P) (2P| 



#2P - ir 2 p/2 - (Eis + H-W/2- (Eis + w L ) 



(27) 



That means that the resonance term involving only the 2P state is replaced by the full Green 
function of the electron, including the off-resonant states and the "decay rate operator" T defined 
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loi rj E 2 p — Eis 




A 2P X S nP X S 2P IS 



(b) 



Figure 3: Same as Fig. 2, but with off-resonant virtual states (nP states with 
The arrow of time is from left to right. 



in Eq. (|24|). The "energy argument" of f = f {E) is to be taken as E = E\$ + lvl ~ in 
Eq. (^), H is a priori the Schrodinger Hamiltonian for the bound electron [we have done all 
calculations in the nonrelativistic (NR) approximation, but the calculations may be generalized 
to the relativistic case]. 

It is now easy to verify by inspection that the amplitudes corresponding to the Feynman diagrams 
in Figs. H (a) and (b) are generated by an expansion of the propagator 1/{H — iT/2— (Eis+ujl)} 
in powers of T. Furthermore, the diagram in Fig. || (b) exactly corresponds to the "second-order 
perturbation" in Eq. (0). So, we conclude that the energy shifts by squared decay rates should 
rather be interpreted as radiative nonresonant corrections to the photon scattering cross section 
than energy shifts of individual atomic levels. Specifically, the "second-order perturbation" in 
Eq. which corresponds to the diagram in Fig. || (b), is merely a nonresonant generalization 
of the resonant diagram Fig. |2] (c). Both diagrams - Fig. Q (c) and Fig. |3| (b) - can be treated 
in a natural way by appropriate replacements within the electron propagator denominator. 
Specifically, we have the replacement H — > H — if/2 for the infinite series involving the diagrams 
in Figs. |H (a) and (b), and the corresponding replacement E^p — > E^p— iT2p/2 for the propagator 
denominator in the resonance approximation [see Figs. [2] (b) and (c) and the right-hand side of 
Eq. (p6|)j. Both added terms (— iT/2 and — ir2p/2) are not real energy shifts. 
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4 Conclusions 



We have investigated a problematic set of two-loop self-energy corrections involving the square of 
the decay rate of the atomic state. These self-energy corrections cannot be interpreted as radia- 
tive energy shifts in the usual sense, although the relevant terms are generated by a fourth-order 
expansion of the Gell-Mann-Low-Sucher theorem (II) in powers of the quantum electro dynamic 
interaction Lagrangian (|3|). As explained in Sec. |, some of the expressions arise naturally if 
we treat the self-energy leading to the decay rate as a first-order perturbation and consider the 
resulting second-order perturbation. The same problematic expressions involving squared decay 
rates result from the alternative formalism for deriving level shifts based on the two-time Green 
function method [19]. At some risk to over-simplification, we can state that the difficulties are 
related to the fact that the current methods inadequately address the question of the preparation 
of the excited atomic states and involve asymptotic states with an infinite lifetime. As explained 
in Sec. ||, some of the problematic two-loop self-energy corrections find a natural interpretation 
as radiative nonresonant corrections to the line-shape for atomic transitions which are of second 
order in the "decay rate operator" T defined in Eq. (|24|). 

For the 2P-level in atomic hydrogen, as shown in Sec. ^ the problematic energy shift by squared 
decay rates is —14.9 Hz [see Eq. (|l~9|)[. For the 3P state of atomic hydrogen, the correction 
amounts to —7.47 Hz [see Eq. (|2l])1. It is perhaps interesting to note that the effect discussed 
here scales as a 2 (Za) 6 mc 2 (for P states and states with higher angular momenta). The same 
order-of-magnitude is characteristic of the shift of the peak of the total photon scattering cross 
section in electric-dipole transitions in atomic hydrogen ]20| — 23 1 . As is evident from Eqs. (16) and 
fll~8|), an order-of-magnitude estimate for the problematic two- loop effect is given by the ratio 
of the square of the decay rate of the atomic state to a typical atomic energy level difference. 
For the 2S state with a natural line width of 1.3 Hz, the squared decay rate therefore does not 
represent an appreciable predictive limit in current and future experiments. The situation is 
different for the 2S-8D transition studied in ]l]] . Parametrically, the "squared-decay" corrections 
are of the order of a 2 (Za) e for the 8D state; however, if we assume a typical l/rz 3 -type scaling of 
the effect (n is the principal quantum number), then we immediately obtain an estimate below 
1 Hz for the problematic effect in the 2S-8D transition in atomic hydrogen. 

There has recently been a dramatic increase in the accuracy to which atomic energy levels can 
be measured experimentally [|| and evaluated theoretically (e.g. |^-[7j ). Tiny nonresonant effects 
that influence the natural line shape of hydrogenic transitions have received considerable atten- 
tion [21-23], and it has been pointed out that nonresonant effects are enhanced in differential 
cross sections as opposed to total cross sections [22,23]. Essentially, the nonresonant effects give 
rise to a process-dependent new line shape, different from the Lorentzian, which has to be fitted 
by a number of parameters. By a suitable fit of the line shape, taking into account properly 
the relevant experimental conditions, it is in principle possible to correct the observed peak of 
the cross section for the nonresonant contributions. The same applies at least in part to the 
problematic two- loop corrections discussed here (see the discussion in Sec. |3]), provided they are 
interpreted properly as radiative corrections to the off-resonance effects on the line shape. The 
problematic two-loop corrections illustrate that it is impossible to separate atomic energy levels 
at the order of a 2 (Za) 6 mc 2 , that is to say to define an energy shift that relates to one and 
only one level. 

In ^4| , it has been stressed that several important problems associated with the normalization of 
electronic states and with the infrared catastrophe can be avoided if we consider the electron in 
a free state as obtained by ionization of a bound state (for example, by ionization of the atomic 
ground state which does not decay at all and is therefore the only "true" asymptotic state). The 
considerations presented in [24] imply that infrared divergences of free-electron quantum electro- 
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dynamics originate because one "disregards the methods of obtaining and detecting these (free) 
states" . It has been shown that in special cases, the spontaneous emission from a certain "well" 
constructed state is exactly exponential p4j. We have mentioned above that the interpretation 
of the energy shift by squared decay rates is problematic because it is calculated as a matrix 
element evaluated on the excited atomic state. Within the formalism introduced in [24|, we can 
argue that this procedure of evaluating the energy shift of the excited atomic state according to 
the current formalism is inconsistent because it disregards the fact that according to [24], the 
only "proper" way of defining the excited state is to view this state as having been obtained by 
excitation from the ground state. 
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